Structural response to 0*-0' and magnetic transitions in orthorhombic perovskites 
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We present a temperature dependent single crystal x-ray diffraction study of twinned orthorhombic 
perovskites Lai-rcCa^MnOa, for x — 0.16 and x = 0.25. These data show the evolution of the crystal 
structure from the ferromagnetic insulating state to the ferromagnetic metallic state. The data are 
modelled in space group Prima with twin relations based on a distribution of the b axis over three 
perpendicular cubic axes. The twin model allows full structure determination in the presence of up 
to six twin fractions using the single crystal x-ray diffraction data. 

PACS numbers: 61.72.Mm 
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I. INTRODUCTION 

The manganites have generated considerable attention 
because of the colossal magnetoresistance effect. The role 
of magnetic order has been widely discussed since double 
exchange, inducing the ferromagnetic order, is required 
to generate a metallic ground state. The role of orbital 
order is much less understood. While local Jahn- Teller 
(JT) distortions are crucial in explaining the localisation 
of the charge carriers in the paramagnetic state, many of 
the structural features of the long range Jahn- Teller (JT) 
ordering has not been well studied and warrant closer in- 
vestigation. Even for undoped systems the JT ordering 
is not well understood as the electronic degrees of free- 
dom have a strong interaction with strain.tH The orbital 
ordering in perovskites with degenerate e g electrons can 
be easily measured, whereas for degenerate t2 9 electrons 
the Jahn- Teller distortions are much smaller. In this pa- 
per we focus on the effects of doping on the JT order- 
ing and the role of the lanthanide site shift. In order 
to minimise the effect of "extraneous" strain, we have 
studied high quality single crystals. A method to obtain 
the crystal structure of twinned orthorhombic samples 
is presented using single crystal x-ray diffraction (SXD). 
Furthermore, we will show how the generally observed 
Pnma symmetry can accommodate a 3D rotation of the 
octahedra, known as the GdFeC>3 rotation, concurrently 
with a Jahn- Teller ordering. 



A. General 

Generally, three basic states for manganites are distin- 
guished, viz. the ferromagnetic metal, the paramagnetic 
polaronic liquid and the orbital and/or charge ordered 
antiferromagnetic insulating state. These states can be 
identified by different structural features. Metallicity in- 
duces charge delocalisation and is associated with equal 



Mn— O bond lengths. The charge and/or orbital ordered 
phases give rise to, for instance, lattice doubling. The po- 
laronic liquid is locally characterised by small polarons. 



B. diffraction 

Typically crystal structure information is generated by 
neutron powder diffraction (NPD). NPD gives direct in- 
formation about the interplanar distances and is very ac- 
curate in determining the changes of the lattice parame- 
ters with temperature and pressure. Using a full pattern 
Rictveld Analysis the atomic positions can be determined 
with an accuracy of 10 -4 -10 -3 . NPD is very powerful for 
oxides because of the large nuclear scattering factor for 
neutrons of oxygen, compared with the low scattering 
factor for X-rays. However NPD is not very suitable to 
determine the symmetry of the crystal structure. 

The method that we have used is SXD. SXD is 
a very powerful experimental application to determine 
crystal structures of many materials. Whereas powder 
diffraction projects the reciprocal space to one dimen- 
sion (d spacing, or diffraction angle 28), SXD retains the 
full 3-dimensional information of the reciprocal lattice. 
The relative co-ordinates of the atoms in the unit cell 
can be determined from the intensity distribution with 
an accuracy better than by NPD. 

However, for many perovskites SXD is not used be- 
cause of the unconventional twinning of the crystals. The 
structure of orthorhombic perovskites is a deformation of 
the structure of a simple cubic perovskite. The deviations 
from cubic symmetry are relatively small. Therefore, do- 
mains with specific orientations of the unit cell can grow 
coherently on each other, a phenomenon known as twin- 
ning. Twinning of the perovskites is complex as the dis- 
tortions from the cubic symmetry allow not only simple 
twin axes and twin planes, but also 3D twin relations. 
Then the observed intensities of the diffraction spots in 
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reciprocal space are the weighted sums of the same lattice 
in different orientations, with a weighting factor propor- 
tional to the volume of the corresponding twin fraction. 
Detwinning is the process of finding the original intensity 
distribution, allowing a full structure refinement. 

A method is presented to analyse the 3D twinned crys- 
tals, and to apply a full structure determination. Be- 
sides it is shown that one does not need a larger data 
set than for untwinned crystals, which is unusual. Using 
SXD we observe a 2a p x 2a p x 2a p unit cell, with a p the 
lattice parameter of the simple cubic perovskite. It is 
commonly accepted that the unit cell of Lai-^Ca^MnOa 
is \/2a p x 2a p x \/2a p . The twinning is unique and it 
involves a distribution of the b axis over three perpen- 
dicular cubic axes.B This leads to a large fraction, about 
25%, of coinciding reflections and an only twice as large 
unit cell. The observed unit cell originates from a three- 
dimensional type of twinning that is not restricted to 
manganites. The model is most likely of general applica- 
tion for a large variety of perovskite Pnma crystals. 

II. CRYSTAL STRUCTURE 

Most structure research of perovskites focuses on the 
Mn-0 distances and the Mn— O— Mn angles, as these 
parameters determine the competition between the su- 
perexchange and double exchange interactions. The ap- 
proach here stresses the importance of a complete struc- 
ture refinement, including the La-position and the rota- 
tion/distortion of the MnC>6 octahedra. The main de- 
formations that determine the deviation from the cubic 
structure are the GdFeOs rotation and the Jahn- Teller 
distortion. 

The basic building block of the perovskite structure is 
a ~ 3.9 A cubic unit cell with Mn in the centre and O 
at the face centres. The oxygen ions co-ordinate the Mn 
ions to form MnOg octahedra. The A atoms are located 
at the corners of the cube. The undistortecL'parent' cubic 
structure rarely exists, for instance SrTiOad, but distorts 
to orthorhombic or rhombohedral symmetry. 

A. GdFeOs rotation 

Due to the small radius of the A site ion, with re- 
spect to its surrounding cage, the MnC>6 octahedra tilt 
and buckle to accommodate the Janthanide ion. This 
is known as the GdFeOs rotationa and yields the space 
group Pnma. The cubic state allows one unique oxy- 
gen position. Due to the GdFeC>3 rotation we need two 
non-equivalent oxygen positions to describe the structure 
in Pnma symmetry. 02 is the in-plane oxygen, on a 
general position, (x,y,z). Two opposite Mn— 02 bonds 
have the same length, but the perpendicular bonds do 
not need to be equal. 01 is the apical oxygen, located on 
a fourfold (x, j, z) position on the mirror plane. Mn— 01 
bonds are always of the same length. Both in the undis- 



torted and the distorted perovskite, the O— Mn— O bond 
angles are 180° (or near 90°), but due to the buckling 
Mn— O— Mn bond angles are significantly less than 180°. 
A pure GdFeOs rotation can be obtained within Pnma 
with equal Mn-0 bond lengths, as is seen in Fig. m. 




FIG. 1: Sketch of the GdFeOs rotation in the ac plane, obey- 
ing Pnma symmetry. Mn and O ions are represented by large 
gray and small open circles, respectively. Open arrows indi- 
cate the movement associated with the GdFeOs rotation. The 
full line indicates the unit cell. 

The GdFeOs rotation is strongly related to the toler- 
ance factor t, which is defined as: 

t= <^ + /3 + )+r 03 - > (1) 
V2((r Mni+ /3+) +r Q 2-) 

where r x is the ionic radius of element x and (r) de- 
notes the average radius. For an ideal perovskite struc- 
ture the ratio between the radii of the A site ion and 
the transition metal ion is such that the tolerance fac- 
tor is equal to one. In the La-Ca system the tolerance 
factor varies from 0.903 to 0.943 going from LaMnOs to 
CaMn03. For the A site, a co-ordination number of nine 
is assumed, values have been taken from Shannon and 
Prewittfl. 

The rotation also depends on the temperature. At 
high temperature all perovskite structures are (near) cu- 
bic, due to the thermal motion of the ions. These mo- 
tions diminish with decreasing temperature, thereby de- 
creasing the 'effective' radius of the ions. As a result, 
the lanthanides prefer a tighter co-ordination and this is 
achieved by an increase in the tilting and buckling of the 
octahedra. 

The rotation is also sensitive to a ferromagnetic order- 
ing. A ferromagnetic ordering is usually accompanied by 
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an increase in the molar volume. For the Prima per- 
ovskites, lessening the rotation will yield an increase in 
molar volume. Therefore, ferromagnetic ordering will be 
characterised by a step in the rotation parameter, with 
the smaller rotation in the ferromagnetic phase. 



B. Jahn-Teller distortion 

The Jahn-Teller effect originates from the degenerate 
Mn 3+ d 4 ion in an octahedral crystal field. Two pos- 
sible distortions are associated with the Jahn-Teller ef- 
fect. Q2 is an orthorhombic distortion, with the in-plane 
bonds differentiating in a long and a short one. Q3 is 
the tetragonal distortion with the in-plane bond lengths 
shortening and the out-of-plane bonds extending, or vice 
versanU. The main result of the distortions Q2 and Q3 
is that the Mn— O distances become different and the de- 
generacy of the t2g and e g levels is lifted. In Fig. ||, a Q2 
distortion is shown. 




FIG. 2: Sketch of the Q2 JT distortion in the ac plane, obey- 
ing Prima symmetry. The arrows indicate the movement as- 
sociated with the JT distortion. The full line indicates the 
unit cell. 



For LaMnC>3, the orbital splitting is such that z 2 -like 
orbitals are occjuiaied, alternately oriented along the x + z 
and x — z axistnJ, yielding the Q2 distortion. The 02 
fractional co-ordinate and this Q2 distortion are thus 
intimately related. For a Q3 distortion, only the ratio 
between the b lattice parameter and the a and c lattice 
parameters has to change. 



C. Glazer's view on the octahedra 

There are three rotational degrees of freedom for a 
rigid Mn06 octahedron. Glazer identified all possible 
sequences of these rotations for the perovskite systcmEJ. 
Two of these sequences create the symmetry elements 
that make uprtiic standard perovskites. In Glazer's nota- 
tion: a~6 + a - liil yields orthorhombic Prima and a~a~a~ 
will result in rhombohedral R3c. As Glazer worked in a 
pseudocubic 2a p x 2a p x 2a p system we have to transform 
the rotation a~b + a~~ to the orthorhombic unit cell, which 
gives a~b + c°. This means that the rotation around the 
c axis, r) z , is zero. Any movement of an oxwen position 
will also affect the position of the La aton£3. 

By rotating the MnOg octahedra only around the 
x axis, two oxygen ions are raised from their cubic posi- 
tion and two are lowered as is sketched in Fig. ||. Here 
we neglect the influence of a rotation r\ y on these co- 
ordinates. This is in agreement with the observation that 
the four positions of 02 near the Mn at (0, 0, 0) have the 
same distance to the y — plane, albeit two are positive 
and two are negative. As a result of r\ z =0 the fractional 
x co-ordinate of 01 is zero. In practice small deviations 
are found, indicative of the non-rigid behaviour of the 
MnOe octahedra. Typically in AMnOa these two rota- 
tions and additionally the Q2 JT distortion are observed. 




FIG. 3: Rotation, rj x , of the MnC>6 octahedron around the 
a axis. The arrows indicate the motion of the 02 ions. 

Another effect of the GdFe03 rotation is that the num- 
ber of formula units per unit cell is enlarged. The unit 
cell is doubled in the b direction, with respect to the orig- 
inal cubic cell. The ac plane is also doubled, resulting in 
b w 2a p , and a w c w \[2a p . Despite the rotations and 
distortion, the structure remains in origin cubic. The 
shifts of the atoms in the unit cell are small. Further- 
more, Mn remains octahedrally surrounded although the 
Mn— O distances tend to differ and some O— Mn— O an- 
gles may no longer be perfectly 90°. The 180° O-Mn-0 
bond angles are intrinsic to the inversion centre at the Mn 
position. 
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All oxygen positions are related by the symmetry op- 
erations of the space group Pnma. That means that the 
in-plane co-ordinates of one ion [xo2 , zoi ) determine all 
other in-plane positions. However, any linear combina- 
tion of these two parameters and its orthogonal comple- 
mentary combination will be as good. Normally we prefer 
X02 and Z02 , since we refine the atomic positions as xoi 
and z<j2- However, the oxygen shifts (arrows) in Fig. [I], 
indicating tbe rotation, can be expressed as a movement 
along [101]E3. We note that these are perpendicular to 
the oxygen shifts resulting from JT ordering shown in 
Fig. | 

Thus, we can express the position of the oxygen atoms 
relative to their "cubic" position by the 'shift' parame- 
ters X02 + Z02 and X02 — Z02 as sketched iiV-Eig. ^. Here, 
a movement of the 02 ion parallel to [101]E3, which im- 
plies that X02 + Z02 — 5, results in equal bond lengths 
but a Mn-O-Mn angle smaller than 180°. We inter- 
pret this movement, along [101], as the GdFeOs rota- 
tion. Similarly, a shift of the O ion along [lOl], thereby 
fixing X02 — Z02 = 0, results in different in-plane bond 
lengths and therefore indicates a Jahn- Teller distortion. 
The Q2 distortion variable is defined as X02 + Z02 — \, 
thus no distortion yields a Q2 distortion parameter equal 
to zero. The undistorted cubic structure obeys both 
X02 + Z02 — \ = and X02 — Z02 = 0. It is clear that any 
02 position can be expressed uniquely as a superposition 
of a GdFeOa rotation parameter, viz. X02 — Z02, and a 
JT distortion parameter, viz. X02 + ^02 — \- 




FIG. 4: Sketch of the GdFe03 rotation (open arrow) and 
the JT distortion (A-arrow) in the ac plane, obeying Pnma 
symmetry. Mn and O are represented by large gray and small 
open circles, respectively. 



III. STRUCTURE ABO3 VS IONIC RADIUS 

The method of rewriting the oxygen position to two 
variables, which indicate the magnitude of the JT order- 
ing and the magnitude of the GdFeOa rotation, is applied 



to two well-studied series of compounds. Marezio et al. 
have studied the structuij&,of the AFe0 3 single crystals 
with SXD in great detail.E-3 Their data allow us to focus 
on the correlation between the ionic radius, ta and the 
magnitude of the rotation. In the top panel of Fig. ^|, we 
see that with decreasing ta the rotation (open squares) 
increases monotonously. We also observe an increasing 
shift from the cubic position for the A site atom with 
decreasing ta- 

The basic argument favouring the A site shift is that 
the band structure energy can be optimised by decreas- 
ing the A— O distances for the three shortest bonds and 
increasing the A— O bond lengths for the three next- 
shortest bonds. This effect is-demonstrated by the cubic- 
tetragonal transition in WO3EJ. The tungsten ion is octa- 
hedrally surrounded by oxygen ions. Band structure cal- 
culations show that the gap is lowered by a ferroelectric 
displacement of the W ion off-centre. The energy is low- 
ered by the overlap between the empty W 5d xz and §dm 
orbitals with the occupied O 2p x and 2p y , respectively^. 

Goodenough showed the importance of interactions be- 
tween the A ion, d°, and the filled, oxygen 2p orbitals in 
the case of a GdFeOa rotationEZrEEl. Recently, Mizokawa 
et al. reported on the interplay between the GdFeOs ro- 
tation, the orbital ordering and the A site shift in ABO3, 
with B=Mn 3+ (3d 4 ) or V 3+ (3d 2 ). Their theoretical cal- 
culations suggest that the observed orbital ordering in 
LaMnOs is stabilised by bath a large GdFe03 rotation 
and a shift of the A site iont3. Conversely, if a Jahn- Teller 
distortion is present, then the energy will be lowered if it 
is accompanied by a shift of the A site. 

In cubic perovskites, the A ion is in the centre of 
the AO12 polyhedra. The GdFeOs rotation distorts the 
AO12 polyhedron. Additionally, minimising the bonding- 
antibonding energy will decrease the shortest A— O bonds 
even more by a shift of the A ion. Marezio et al.'s data 
show-the correlation between the rotation and the La site 
shiftli 4 !. The relevant atomic parameters, x and z (open 
circles) of the A site, and the rotation parameter (open 
squares) are fully correlated. The bottom panel shows a 
perfect linear relation between the A site shift and the 
rotation parameter of AFe03 (open symbols). Thus the 
shift of the A site atom is, in this system, fully deter- 
mined by the rotation of the FeOe octahedron. This con- 
clusion is supported by inspection of the neutron powder 
diffraction data on AMn03 (closed symbols) by Alonso 
et aLu-3 Their data show that the JT distortion param- 
eter is ~ 01)301(8) at T = 295 K, independent of the A 
ion radiusEj. The observed slope of rotation vs. site shift 
is identical for AMn03 and AFe03 within the accuracy 
of the measurement. 

The A site shift is the result of the covalent interaction 
between the A 3+ ion and the surrounding oxygen ions. 
The data for AFe03 shows that the ratio xa/za = 5. 
Thus the A ions will decrease their distance to the near- 
est three oxygen ions by moving along a vector [1 
5]. As the A co-ordination is directly determined by 
the changes in rotation and distortion, the A site shift 
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FIG. 5: Top panel: The rotation and zj^ versus the ionic ra- 
dius of AFeC>3 and AMnC>3. Both parameters indicate a de- 
viation from the cubic perovskite, and both increase with de- 
creasing rj^. Bottom panel: The rotation versus z^. The pa- 
rameters show almost perfect correlation indica ting their in- 
timate relation. Data have been taken from ref. jFf (AFeOs) 
and ref. |l^] (AMnOa). Open and closed symbols represent 
AMn03 and AFe03, respectively. Squares and circles repre- 
sent the rotation and site shift parameter, respectively. 



strongly depends on both the GdFeC>3 rotation and the 
JT distortion. 

Although thorough structural studies on the mangan- 
ites are published in the literature, see e.g. ref. |l9[ they 
typically are based on NPD data of Lai-^Sr^MnOa. In 
1971, Elemans et al. studied the structure of LaMnOa 
and Lao.95Cao.o5Mn03, using NPD. They found that 



the Q2. distortion is similar for T = 4.2 K and T — 
298 KB The Mn-0 bond lengths for La .75Cao.25Mn0 3 
are reported by Radaelli et al. They show a reduc- 
tiop-in the distortion or bond length disproportion at 
T c .E!l Recently, Huang et al. report the Mn-0 bond 
lengths for Lao.ssCao.isMnOa, Lao.7sCao.25Mn03 and 
Lao.67Cao.33Mn03 using NPD. The scatter and error bars 
make any distinct observation hard. A slight decrease in 
variance of the bond lengths is observed near T c El In a 
Raman study, the MnOe distortion was reported to de- 
crease linearly from 42% at x = to 8% at x = 0.4, using 
x-ray powder diffraction.E3 Dabrowski et al. studied the 
relation between oxygen content, doping concentration 
and crystal structure for Lai-^Ca^MnOs. They report a 
transition at room temperature from the O' state, coher- 
ent JT distortions, to the O* state at x = 0.14.E2lNo sign 
of coherent JT distortions was found for Lao.gCao^MnOs 
and Lao.7Cao.3Mn03 down to T = 15 K. Increasing the 
temperature towards T c , increases the variance in bond 
lengths, as was determined from total correlation func- 
tion experiments. This however measures the local JT 
distortions .cJ 



IV. EXPERIMENTAL PROCEDURES 

The experiments were carried out on single crystals of 
Lai-zCa^MnOs, x = 0.16 and x — 0.25, obtained by the 
floating zone method at the JRCAT institute, Tsukuba, 
Japan. Although all crystals were twinned, small mo- 
saicity and sharp diffraction spots were observed. Fur- 
thermore, sharp magnetic and electronic transitions in- 
dicate the good quality of the crystals. A thin piece 
was cut from the crystal to be used for single crystal 
diffractomctry. The single crystal was mounted on an 
ENRAF-NONIUS CAD4 single crystal diffractometer. 
The temperature of the crystal was controlled by heat- 
ing a constant nitrogen flow. Initial measurements were 
done at 180 K. Temperature dependent measurements 
between 90 K and 300 K were performed on a Bruker 
APEX diffractometer with an adjustable temperature 
set-up. Lattice parameters based on observations by the 
CCD camera are refined using all observed reflections. 
Details of the SXEj— experiments and refinements are 
published elsewhere £j The single crystals exhibit sharp 
magnetic transitions with T c ~ 160 K and 200 K, for 
Lao.84Cao.i6Mn03 and Lao.75Cao.25Mn03respectively. 



V. TWINNING 

A. The twin structure for manganites 

We will focus in this section on the main reason why 
SXD has not been widely used for these perovskites, 
namely twinning. Twinning in doped LaMn03 originates 
from the transition of the highly symmetric cubic parent 
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structure to the orthorhombic symmetry, that accommo- 
dates both the GdFeC>3 rotation and the Jahn- Teller dis- 
tortion. The solution of the twin relations in the crystals 
allows us to study in detail the ordering of these com- 
pounds as influenced by temperature, magnetic state and 
doping concentration. But for now we will focus on the 
twin relations. 

Twinning can complicate the structure determination 
and is often observed in crystals with a reduced sym- 
metry. The orthorhombic perovskite LaMn03 has lower 
symmetry than cubic SrTiC>3. Conventional twin mod- 
els keep one characteristic axis unchanged and form the 
domains by rotation around that axis. This is com- 
monly observed in constrained epitaxial thin films and 
pseudo-two-dimensional crystals like YF^CusCv-a. An- 
other standard twin is inversion twinning, found in non- 
centrosymmetric systems. Due to the inversion twin, 
pscudo-centrosymmetry is obtained and the net p plarif - 
sation in ferroelectric compounds is reduced to zerocSEZl. 
Usually, both types of twinning have the property that 
either all twin domains have reflections that lie on top 
of each other, merohedral twins, or there is a non- 
commensurate set of reflections, e.g. non-merohcdral 
twins. Here a more extensive form of twinning is pro- 
posed, where at a quarter of the observed reflections the 
reflections of different domains coincide and the other 
part of the observed reflections originates from a single 
domain only. Together they give rise to the observed 
metrically cubic lattice. 

Rodrigruez-Carvajal et al. noted that their neutron 
powder spectra of pure LaMnC>3 could be indexed in a 
cubic 2a p x 2a p x 2a Pr -unit cell, above the Jahn- Teller 
transition tcmpcraturcES. They did not observe an or- 
thorhombic splitting of the peaks. Nevertheless, they 
could not refine the structure in a cubic unit cell. But full 
pattern refinement was-only possible in the orthorhom- 
bic space group PnmaEB. Note that in powder diffraction 
data there is no direct method to observe the threefold 
symmetry along the body diagonal. The refinement has 
been done in the conventional Pnma setting. They sug- 
gested that the powder is twinned but could not give 
evidence, as they only studied powders. Combining our 
observation, including the double cubic lattice, the sys- 
tematic extinctions and the non-cubic intensity distribu- 
tion, with the suggestion given by Rodrigruez-Carvajal 
et al. leads to the conclusion that our single crystals are 
twinned. 

Electron diffraction and microscopy are instrumental 
to characterise this twinning because of its high spatial 
resolution. Recently Dechamps et aln identified the twin 
relations of orthorhombic and rhombohedral AMn03. 
Their data are interpreted for Pnma symmetry with co- 
herent domains, with the doubled b axis in three orthog- 
onal orientations. 

Reflections in reciprocal space are experimentally ob- 
served at a regular distance in three orthogonal direc- 
tions, corresponding to a cubic lattice spacing of 7.8 A in 
real space. Although the three axes had equal lengths, 



we could not observe the required threefold symmetry 
axis along (111) for cubic symmetry. This is shown for 
eeo reflections in Fig. ^, E denotes an even value for the 
Miller indices, O stands for odd. Furthermore, study- 
ing the intensity distribution of hkl planes with constant 
h, k or I showed much regularity, and some anomalous 
extinction conditions are observed. An overview of these 
features is given in appendix |a|. A twinned structure con- 
sisting of coherently oriented Pnma domains isjproposed, 
analogous with the results of Dechamps et aln. Due to 
partial overlap this results in a metric cubic system with 
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FIG. 6: The number of reflections versus the orientations OEE, 
EEO and EOE for different intensity ranges. OEE, b' parallel 
to a" , has mostly very low intensity. Reflections EOE, b' par- 
allel to 6", have a broad range of intensities. In cubic sym- 
metry, these distributions should be identical for absorption 
corrected data, which we considered. 



B. Twin model 

The transformation of cubic to orthorhombic symme- 
try requires a designation of a, b and c with respect to 
the degenerate cubic axes. There are three possibilities to 
position the doubled b axis along the three original cubic 
axes. Thus we propose that the three fractions' b axes are 
oriented along the three original axes of the cubic unit 
cell, as sketched in Fig. 0. This twin model consists of a 
Pnma unit cell, transformed by rotation about the 'cu- 
bic' [111] axis. We still have the freedom to choose the 
a and c axes, perpendicular to the b axis, and rotated 
45° with respect to the cubic axes. Therefore, this model 
yields six different orientations of the orthorhombic unit 
cell. As the differences between a/y/2, 6/2 and are 
small, we observe the reciprocal superposition of the six 
orientations as metrically cubic. 

Using the proposed model, the different contributions 
of the twin fractions to the total integrated peak inten- 
sity are taken into account, according to the following 
procedure. We considered for each observed reflection 
the corresponding orthorhombic reflections for the twin 
fractions (if applicable) and used an identical unit cell 
for each twin. The refined model consisted of the regu- 
lar parameters in single crystal diffractometer. This way, 
both the atomic positions of the asymmetric unit and the 



b"=b 





FIG. 7: Two of three possible orientations of the Pnma unit 
cell in the observed 2ap x 2ap x 2ap (a" x b" x c") unit cell, 
bi , &2 an d &3 denote the orientations of the doubled axes along 
the original cubic axes, ap. Note that the choice of a and c 
is still open, after choosing an orientation for the b' axis. 



volume ratio of the twin domains are refined simultane- 
ously. 

Several crystals were measured and refined with this 
detwinning model. The refinements showed that for ev- 
ery crystal the distribution of the volume over the twin 
fractions is different. This is another indication that our 
crystals are twinned. Other structural deformations than 
twinning would lead to a constant, not sample depen- 
dent, effect on the structure factors. Attempts to refine 
will then result in the wrong model, giving the same value 
for the volume fractions, independent of sample. We can 
also conclude that the size of the twin domains must be 
smaller than the magnitude of the measured crystals, i.e. 
0.1-0.2 mm. Larger twin domains would give rise to crys- 
tals of one single domain and even smaller domains are 
more likely to produce a constant spreading of the do- 
mains. Non-regular, though constant, volume fraction of 
the twin domains could also signal a preferential growth 
direction of the crystal. 



VI. RESULTS 
A. Lattice parameters 

Reflections in reciprocal space are observed on a dou- 
ble cubic lattice. From these lattice spacings, only the 
'average' lattice parameters can be calculated. In Fig. || 
the lattice parameter a is plotted against temperature. 
For Lao.84Cao.i6Mn03 the volume fractions of two of the 
three b axis orientations were very small, therefore, ef- 
fectively a normal pattern is observed, within only the 
a-c/c-a twin. For this sample both lattice parameter a 
and b are plotted. Lattice parameter c behaved identical 
to lattice parameter a. 
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FIG. 8: The temperature dependence of the observed lattice 
parameters a. For Lao.84Cao.i6Mn03 lattice parameter b is 
plotted on the right axis. The error bars are of the same size 
as the symbols. Note that the error bars indicate an under 
limit to the uncertainty. 

The lattice parameter for Lao.75Cao.25Mn03 shows a 
continuous increase with temperature. No anomaly near 
T c can be observed. Lattice parameters a and b for 
Lao.75Cao.25Mn03 respectively show a sharp down turn 
and up turn above T ~ 250 K. 



B. Rotation 

In Fig. ^ the temperature dependence of the rota- 
tion is plotted versus the temperature. The rotation- 
temperature curves depend strongly on the Ca concen- 
tration. The curve of the x — 0.25 sample shows a sharp 
steplike increase in the rotation at T c . At x — 0.16, the 
scatter in the data is larger, though a smaller step near 
T c can still be observed. Above T c , the rotation decreases 
with a similar slope as the curve for x = 0.25. The ab- 
solute value of the rotation is about 10% larger for the 
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lower doped sample. 
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FIG. 9: The temperature dependence of the rotation param- 
eter. The broken lines are guides to the eye. 



C. JT distortion parameter 

Fig. |Io| shows that the sample with x = 0.16 has a con- 
stant value for the JT distortion parameter at low tem- 
perature. Above T = 260 K the JT distortion parameter 
decreases rapidly with temperature, which signals the JT 
orbital ordering temperature. 
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FIG. 10: Jahn- Teller parameter against temperature. The 
broken lines are guides to the eye. 



position (0,4,5). W e note that the error bars of the 
A position, 0.0001, and scatter are much smaller than the 
error bars on the JT distortion parameter and rotation 
parameter, 0.0005. These latter are derived from the 
oxygen positions, which have a larger uncertainty due to 
their smaller electron density. 

The ratio of xa and za is constant, with xa roughly 5 
times larger than za for all temperatures. 
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FIG. 11: The temperature dependence of the A site position 
of Lao.84Cao.i6Mn03 (top) and Lao.75Cao.2sMn03 (bottom). 
The ratio between xa and za is 5.6 and 4.8. The xa axes 
have the same scale to allow comparison. The corresponding 
za axes have been scaled with the average ratios 5.6 and 4.8. 



VII. DISCUSSION 



A. Lattice Parameters 



The JT distortion parameter for the x — 0.25 sample 
is independent on temperature within the uncertainty of 
the experiment. The value of 0.002 results from the fact 
that the 02 position is not constrained by symmetry. 



D. La site shift 

In Fig. O, we show the temperature dependence of 
the A position, (xa, \ ,\~ z a), with respect to the ideal 



The occurrence of twins in our single crystals effec- 
tively averages the lattice parameters a, b and c, since 
reflections with different hkl, originating from different 
twin fractions, coincide in one reflection spot. However, 
for the x — 0.16 sample, only twin fractions with par- 
allel b axes are observed, which allows a separate deter- 
mination of a/c and b. The lattice parameter c is still 
"averaged" with a. 

The lattice parameters, or volume of the unit cell, are 
influenced by several effects. The most obvious in a tern- 
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perature dependent study is the thermal expansion. We 
clearly observe the continuous increase in lattice param- 
eter a with temperature for x — 0.25. The GdFeC>3 rota- 
tion is another non-volume conserving distortion of the 
lattice. The main reason the MnC>6 octahedra are rotated 
is the ionic radius of the A site atom. We expect a change 
of volume at T c , but we could not observe an anomaly 
in the lattice parameters at T = 200 K (x — 0.25) or at 
T = 160 K (x = 0.16). 

An orbital ordering can have an effect on the lattice 
parameters, even though it is volume conserving. The 
effect is different for the Q2 and the Q3 type distortion. 
The coherent Q2 distortion does not change the lattice 
parameters in first approximation as the long and short 
Mn-0 bonds are oriented alternately, both in the [101] 
and the [10T] direction. However, the Q3 distortion has 
a significant effect on the lattice parameters. Because it 
has an opposite effect on the Mn-Ol bonds, parallel to b, 
and the Mn-02 bonds, in the ac plane. 

Note that there is a large scatter in Fig. ||. Combined 
with the averaging effect of the twin domains, the re- 
sponses of the crystal structure to the orbital ordering 
or the ferromagnetic transitions are not very clear. This 
stresses the importance of an alternative way to study the 
relation between the structure and the electronic proper- 
ties. 



B. Rotation 



C. JT distortion parameter 

In Fig. [l(| we observe that the Q2 distortion for the 
x = 0.25 sample is constant with temperature, but not 
equal to zero as expected for the ferromagnetic metal- 
lic phase. The value, 0.002(1), is in good agreement 
with values calculated from structure determinations of 
AFe0 3 , e.g. 0.0036 for LuFe0 3 .y Fe 3d 5 in the high spin 
state, has a completely filled subband, and is therefore 
not sensitive to the JT effect. However, other effects, for 
example related to the hybridisation associated with the 
A site shift, and the degrees of freedom of the 02 posi- 
tion in Prima symmetry allow the Q2 distortion to be 
non-zero. 

The reduction in Q2 parameter, as shown by the curve 
for x = 0.16, clearly coincides with the JT ordering tem- 
perature. As single crystal diffraction is sensitive to long 
distance correlations, this either means that the orbitals 
becomes disordered or that there is an overall reduc- 
tion in the magnitude of the local Q2 distortion. Local 
probes, like pair distribution functions, show that even 
in the paramagnetic phase the local distortions do not 
disappear .Ell Our observation of a decrease in Q2 distor- 
tion is therefore caused by a decrease in the coherence of 
the orbital order. 

The data for LaMnOa show that the JT distortion pa- 
rameter is about 0.0320(10) at room temperature and 
that the JT ordering reduces to 0.0060 at T fh?Q8 K, 
where the orbital order is known to disappear.ESEII The 
low temperature JT distortion parameter and the JT 
transition temperature strongly depend on the doping 
level, x. With increasing doping both variables decrease. 



The rotation parameter is clearly sensitive to the ferro- 
magnetic ordering in the x = 0.25 samples at T c — 200 K. 
Likewise a smaller transition is shown in the x — 0.16 
curve in Fig. |9| The decrease in rotation parameter indi- 
cates a straightening of the Mn— O— Mn bond angles as 
shown in Fig. [I] and the text. If a constant volume of the 
MnOg octahedra is assumed, the decrease in rotation, 
associated with the paramagnetic to ferromagnetic tran- 
sition, yields an increase in total volume. An increase in 
molar volume, or an increase in the lattice parameters is 
a common feature of ferromagnetic ordering, as observed 
for instance in pure nickelo This agrees very well with 
our observation of the decreasing rotation parameter and 
consequently increasing volume of the unit cell. 

The rotation for LaMn03 is 0.0810(2Lat room temper- 
ature and decreases with temperature.ES Elemans et al. 
reported no significant change in rotation between room 
temperature and 4.2 K, but their resolution is smaller. 
However, by doping with 5% Ca a decrease in the ro- 
tation parameter to 0.077(1) is observed.EI CaMn03 is 
(almost) cubic and consequently has no rotation. There- 
fore we expect a continuous decrease in the rotation with 
increasing doping. This agrees well with our data, which 
show a decrease of the absolute value of the rotation with 
5-10% between x = 0.16 and x = 0.25. 



D. La site shift 

The ratio of x a and za corresponds to a shift of the A 
atom along [1 5] _ Mizokawa et aZ.assumed the A site 
shift along [1 7] .El Their assumed [1 7] is in good 
agreement with our experimental observation of [1 5]. 

In fig. [H], changes in the gradient are observed near 
T c and the orbital ordering temperature, corresponding 
with the step in the resistance curve. This supports the 
assumption that the A site shift is sensible to changes 
in the oxygen environment, by changes in the GdFe03 
rotation and the JT distortion. 

The A site shift, expressed by xa, in LaMn03r 
0.0490(2) at room temperature and 0.0217(3) at 798 Kb 
The latter value is comparable with the observed values 
near 300 K for x = 0.25. 

If we compare these curves with the curves in Fig. ^ 
and Fig. 10 it is obvious that the curve of rotation vs. 



temperature of the x = 0.25 sample is mimicked by the 
A site shift vs. temperature curve. Similarly, the A site 
shift vs. temperature curve of the x — 0.16 sample shows 
the same rapid fall above T — 260 K as the curve of the 
JT distortion parameter vs. temperature. But a small 
change in slope is observed at T c . Evidently the A site 
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shift is a good probe to pinpoint transitions which involve 
changes in the oxygen co-ordination of the Mn ions. 



VIII. CONCLUSION 

We conclude that full structure determination using 
single crystal x-ray diffraction provides accurate infor- 
mation regarding the magnetic and orbital ordering in 
manganite perovskites. We showed that the twinning 
of these single crystals is based on a distribution of the 
doubled b axis over three perpendicular cubic perovskite 
axes. The shift of the A site atoms is very sensitive to 
changes in the oxygen environment of the A site. There- 
fore, the A site shift is an accurate probe for magnetic 
transitions and orbital ordering, since these involve rota- 
tions and deformations of the MnOg octahedra network. 
This warrants structure investigations close to the fer- 
romagnetic insulator to ferromagnetic metal transition 
near x — 0.2. 
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APPENDIX A: DATA ANALYSIS 

Although cubic symmetry could not be observed, 
some regularities and anomalous extinction conditions 
are present in the intensity distribution over hkl space. 
In Fig. [l^ the observed intensity distribution of the re- 
flections is plotted. The reflections are subdivided in four 
groups, on the basis of the number of odd indices in the 
double cubic setting. Here the effect of the twin model 
on the visibility of the standard Prima reflection con- 
ditions is described and a detailed analysis of the inten- 
sity distribution of a model crystal with Pnma symmetry 
and of the measured reflections in double cubic setting 
is given. The first step in space group determination 
is to search for systematic absences in the list of reflec- 
tions. We observe no intensity for the reflections shown 
in Table |. Note that the indices are in double cubic 
setting and therefore allow cyclic permutation, i.e. hOO, 
OkO and 001 are all represented by hOO. These reflection 
conditions only allow the space groups P2i3 and P4232. 
Attempts to determine the structure using these space 
groups were unsuccessful. But, these space groups only 
have hOO : h = 2n as reflection condition. This suggests 
that we have more information than can be attributed to 
these space groups. Furthermore, the observed and, for 
cubic systems, anomalous reflection conditions suggested 
that we might be looking at a twinned crystal. In the 



TABLE I: Extinct reflections as observed in the 2ap x 2ap x 
2ap data set. 

Reflection Extinction condition Observed reflections 



h" h^2n 

ti'hFb h^2n 

h" k" h" h=£2n,k = 2n 

h" k"hP h=£ 2n, k = 2n 



e00 3 2 

EEOp 

EOE i2 
EOE 



next paragraph, we present the transformation of the re- 
flection conditions for orthorhombic Pnma according to 
the presented twin model. 




I? 



1 



Intensity group 



FIG. 12: The intensity distribution of reflections. Low 
intensities,/ < 200 are mostly EOO, while high intensities, 
I > 2500, are of the types eee and eeo, e denotes an even 
value for the Miller indices, O stands for odd. 

The following equations are used to transform the sim- 
ple cubic Miller indices hkl to the standard Pnma setting 
and to transform these to the double cubic setting. 



h'k'l' = h + l, 2k, h-l 
h"k"l" = h' + I', 2k', ti - V 



(Al) 
(A2) 



The reflection conditions for Pnma are shown in Ta- 
ble O. Using Eq. A2 the Pnma reflection condition 
k : k' = 2n is transformed to double cubic as 
k" : k" — 2n and with cyclic permutation yielding 
h" : h" = 2n and I" : I" = 2n. This corre- 
sponds to the reflection condition h" : h" = 2n in 
Table Q. However, Pnma reflections of the type h! h! 
are also transformed to h" 0, with h" = 2h! . Note that 
the orthorhombic to double cubic transformation gives 
Q k' Q = h" 0, with b' || a", and the reflection con- 
dition k' = 2n implies that only h" = 2n are present. 
While the h! h! reflections intrinsically give h" with 
h" = 2n, with b' || b" . This implies that h' h' may 
never contribute intensity to reflections h" 0, with h" 
odd. Conversely, if we do not observe intensity on reflec- 
tions h" : h" 2n, then all of the constituting Pnma 
reflections, at most six, should be absent. 
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TABLE II: Reflection conditions for space group Pnraa. 



Okl 


k+l 


= 2n 


OkO 


k = 


2n 


00/ 


I = 


2n 


hkO 


h = 


2n 


hOO 


h = 


2n 



To transform and understand the other reflection con- 
ditions of Pnma we have to take into account that an 
orthorhombic, and also a cubic, unit cell implies that hkl 
is equivalent upon sign reversal of each of the Miller in- 
dices. The Pnma reflection condition hi : h' = 2n 
becomes h" h" : h" = 2n and similarly 0/' transforms 
to I" I". But as we just stated, h" h" is identical to 
h" h" . If we combine reflections then a reflection con- 
dition is fulfilled if any of the contributing parts has in- 
tensity. To obtain extinction, all contributing reflections 
should be extinct. As /' has no extinction condi- 
tion, i.e. intensity for all /', they will contribute to all 
h" h" I h" h" reflections, and the Pnma reflection 
condition h! : hi — 2n is masked. 

The overlap of different reflections also occurs for 
k' h! and h! k' 0. These transform to h" k" F 7 , and 
h" k" hi' . Therefore, we cannot disentangle k' hi and 
hi k' in the measured double cubic cell. The Pnma 
reflection conditions for k' V and hi k' are k'+V = 2n 
and hi = 2n, respectively, h" k" h" and h" k" h" are 
found to be extinct for h" ^ 2n and k" = 2n in the 
double cubic setting. To have these reflection extinct, all 
contributing Pnma reflection conditions must be extinct. 
This implies that both hi — 2n and k'+h! = 2n, which 
are the reflection conditions for the contributing k' hi 
and h! k' 0, may not be fulfilled. Thus, taking into ac- 
count that h" = h'+O and k" = k', it is easily derived 
that h" 7^ 2n and this yields k" — 2n. This means OEO 
reflections are extinct. 

Now we consider reflections of the type h" k" h" that 
are not extinct. Our twin model yields six fractions that 
can contribute to the considered reflections. 

1. h" — 2n and k" — 2n. Reflections eee have con- 
tributions from the four fractions with the doubled 
b' axis parallel to either a" or c", this is also true 
for reflections OOO. In addition, we have contribu- 
tions from the two fractions with b' parallel to b", 
i.e. from h! kl and k' hi. 

2. h" 7^ 2n and k" ^ 2n. Reflections OOO, in addi- 
tion to the four fractions with the doubled b' axis 
parallel to either a" or c", have one additional con- 
tribution from k' h! . In the special case that 
h" = k", for both eee and OOO, then two of the 
four fractions, parallel to a" and c" contribute. 

3. h" — 2n and k" ^ 2n. Reflections eoe, only have 
a contribution from the hi k 1 reflection, with b' 
parallel to the observed b" axis. 



We have shown that the intensity of observed h" k" h" 
reflections depends on odd or even. We conclude that 
h" k" h" reflections are extinct for OEO, while the other 
reflections have contributions from up to six twin frac- 
tions. 

We can also learn something about the twinned ori- 
gin by studying the intensity distribution. Our twin 
model provides a natural explanation for the hierarchy 
in intensities of the reflections. First, we consider the 
known transition from the single cubic unit cell to the or- 
thorhombic one to elucidate the patterns that are inher- 
ent to perovskites. Then we proceed to show how these 
patterns are convoluted in the double cubic, twinned 
model. 

The intensity distribution of a normal Pnma per- 
ovskite is shown in Fig. 13, In conventional orthorhombic 
Pnma the reflections hl~K V with h'+V = 2n are stronger 
than h'+V 7^ 2n. From these, the k = 2n reflections are 
stronger than the k ^ 2n ones. This can be attributed 
to their origin in the 3.9 A cubic structure. The hi k' V 
reflections in orthorhombic setting or igin ate from hkl in 
the single cubic cell according to Eq. (Al). Therefore, hi 
and V will always be both even or both odd if the reflec- 
tion hi k' V is related to an integer single cubic crystal 
plane. In other words if hl-V ^ 2n then the reflections 
originate from a crystal plane in the single cubic set- 
ting with a non-integer Miller index, i.e. a superlattice 
reflection and their intensities are typically weak. The 
same argument can be applied to k. As kl = 2k, reflec- 
tions with k' — 2n originate from a regular crystal plane 
and reflections with k' ^ 2n are superlattice reflections. 
These reflections with kl ^ 2n however appear to have 
somewhat more intensity than the hl-V ^ 2n reflections. 

We observe the following patterns in the intensity dis- 
tribution of the double cubic cell. First, reflections eee 
are generally the strongest, eooE3 usually absent, or very 
weak. The following hierarchy can be made if we dis- 
regard the cubic symmetry, eoe reflections were much 
stronger than eeo, which are again stronger than OEE. 
Likewise, OEO reflections were, although weak, stronger 
than OOE and again stronger than EOO. 

In conventional Pnma structure we found that if 
h'+V 7^ 2n then reflections with k' 7^ 2n are stronger 
than those with k' — 2n. In our even/odd nota- 
tion, for the Pnma structure this yields (eee=OEO) >> 
(ooo=eoe) >> (ooe=eoo) > (oee=eeo). So there 
are roughly four intensity groups. We can correlate these 
four groups with the four types of reflections we measured 
in the super cubic setting, with the following intensity hi- 
erarchy: eee>>eeo>>ooo>eoo. 

Here we consider the transformation from the or- 
thorhombic to the double cubic setting for the in-plane 
indices, hi and /'. Double cubic indices are calculated 
with h'+V and hl-V , therefore hi' and I" will be both even 
or both odd. Reflections of the type E?0 or 0?E will have 
no intensity of the twin with b' parallel to b" , as they are 
not originating from an 'integer' crystal plane. In con- 
trast with the single cubic to orthorhombic Pnma transi- 
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FIG. 13: The calculated intensities for a conventional per- 
ovskite, subdivided with respect to odd and even of the Miller 
indices, k — 2n reflections are strongest, if b! + I' = 2n. Re- 
flections with h'+l' 7^ 2n have roughly 100 times less intensity. 



This suggested that one twin had 70% of the volume, the 
others 20% and 10%. 

These reflections are sorted to find the correspond- 
ing orientations. The sort parameter is the observed 
double cubic axis that corresponded with the odd in- 
dices in OEE, and the even indices in EOO. Only cy- 
cled reflections that occurred three times with measur- 
able intensity were taking into account. We found that 
I a ii : lyi : I c n = 5 : 80 : 15. The constraint that all three 
intensities have positive values, ignores the weakest re- 
flections. This method also ignores Bijvoet pairs, as we 
only measured —20 < I" < 2. In principal EOO and EOO 
should have identical intensity. If we indeed had a cubic 
system then these variations should be zero within stan- 
dard deviation. The distribution of the intensities, with 
respect to the different axes, strongly suggests that the 
crystal is twinned. 



tion, where orthorhombic superlattice reflections are al- 
lowed, the orthorhombic to double cubic transformation 
is just a new choice of reference system, that may not 
and cannot generate new reflections. But we have a mix- 
ing of different orientations of the b' axis along the three 
cubic axes a", b" and c", which do generate reflections 
where none are observed without the twinning. We can 
cycle the indices of these reflections so that we get either 
e?e or 0?0. From the corresponding twin, with b' paral- 
lel to either a" or c" , we do have intensity on these e?0 
and 0?E reflections. EOE and OEO will only have con- 
tributions from this particular setting with b' \\ b" . eee 
and OOO can be cycled and in general have contributions 
from all six orientations. 



a. intensity 

The small deviations from the cubic symmetry are re- 
flected in the intensities of the reflections. To illustrate 
this point, we transform an arbitrary, su bcel l, reflection 
hkl to the orthorhombic setting by Eq. (Al). Thus or- 
thorhombic reflections with k' even and h! and V both 
even or both odd stem from planes that already existed 
in the original cubic unit cell. We observe in all diffrac- 
tion patterns of orthorhombic perovskites that reflections 
that originate from cubic crystal planes generally have 
higher intensities than those have that do not. 

The intensity of a particular twin fraction is propor- 
tional to the volume of that twin fraction. Using the 
uniqueness of some reflections, i.e. belonging to only one 
twin fraction lattice, we can search for the orientation of 
the largest fraction. All b' axes are oriented along one 
of the measured double cubic axis. It is impossible to 
differentiate the two fractions with parallel b' and per- 
pendicular a 1 and c'. We considered groups of reflections 
with cycled indices, having either one odd, OEE, or one 
even index, EOO. Roughly, the intensities of the reflec- 
tions within these groups occurred with ratio 70:20:10. 



1. Intensity distribution 

Careful analysis of the intensity distribution shows ex- 
tra evidence for the proposed models. Two important 
observations are made. First, the full data set shows 
that reflections hkl of type eee are by far the strongest. 
Secondly, reflections having one odd Miller index, OEE, 
EOE nband eeo, are second strongest in intensity. Of 
these reflections eoe has the highest intensity as seen in 
Fig. §. 

To understand the observed intensity distribution in 
the 2a p x 2a p x 2a p unit cell, we examine the influence 
of the structural deformations on the structure factors. 
Despite the effect of Jahn- Teller distortions, GdFe03 ro- 
tations and twinning the structure remains in origin cu- 
bic. This means that the intensity distribution of the 
peaks with high intensity will always mimic the inten- 
sity distribution of the peaks of the undistorted cubic 
perovskite. An arbitrary reflection, in the small cubic 
subcell, hkl, is tra nsformed to the orthorhombic unit cell 
using equation Al Thus orthorhombic reflections with h' 
and I' both even or both odd and k' even stem from 
the original cubic unit cell. These reflections therefore 
have the highest intensity, see Fig. [ll| Although the fact 
that the h! k' V — eoe reflections do not originate from 
the original cubic planes, they are allowed in the Pnma 
symmetry. In Table |n] the reflection conditions for the 
Pnma space group are given. Extinctions observed in the 
double cubic setting are thoroughly discussed before. Re- 
flections hi k' I' that satisfy the condition h' k' I' = eoe 
have higher intensities than those that do not. 

With these observations, we return to the observed 
double cubic lattice. To double the Pnma unit cell to a 
2a p x 2a p x 2a p unit cell, a reflection h' k' V is transformed 
to h" k" I" using Eq. (|A2|). The transformation yields 
that double cubic reflections with k" — 2n, k" — k', orig- 
inate from single cubic reflections with k = n. Likewise, 
k" = 2n + 1 reflections originate from non-integer sin- 
gle cubic indices k, and will therefore have less intensity. 
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This explains why the measured h"k"l" =eee reflections 
are strongest and eoe are one order of magnitude less 
strong. 

We have seen that ft" = h'+V and I" = h'-V thus 
ft" and I" are always both even or both odd. Therefore 
eoe only has contributions of the fractions with k" cor- 
responding to the orthorhombic b axis. OEE and eeo 
originate from fractions with their b axis parallel to a" 
and e", respectively. The fact that our observed eoe is 
stronger than OEE and eeo, indicates that, by chance, 
the observed b" axis is parallel with the b' axis of the 
largest twin fraction. 

2. Refinement 

Both the Ca content and the volume fractions of the 
twins can be determined from single crystal x-ray diffrac- 
tion in the following way. First a full hemisphere of re- 
ciprocal hkl space is measured. On this very large data 
set, the atomic positions, the Ca content and the vol- 
ume fractions are refined. The results of the refinement 
are compared with a refinement on only one octant of 
these measured reflections. The outcome was equivalent 
within the error bars. We refined again fixing the twin 
fractions and Ca concentration to the values found for 
the first refinement on the large data set, and allowing 
only the atomic co-ordinates and anisotropic displace- 
ment parameters to change. This resulted in the same 
structural model as found with the complete structure 
determination on the largest data set. 

It is concluded that apparently the refinement is in- 
sensitive for the size of the data set. However, it will 
be dependent on the 9 range of the data set. The in- 
sensitivity can be understood if we view the symmetry 
relations and Friedel pairs. Due to the Patterson symme- 
try, mmm, the relation F(hkl) = F(hkl) is valid for all 
structure factors. Friedel pairs are reflections that have 
intrinsically the same structure factor, related thus by 
symmetry and ignoring absorption and anomalous scat- 
tering effects. Although we have six different orienta- 
tions, one octant chosen in a particular orientation will 
still contain data of all fractions. The six orientations can 
be transformed in such a way that every h' k' V is repre- 
sented for all three b directions in one octant of reciprocal 
space. 

APPENDIX B: IMPLEMENTATION IN 
REFINEMENT 

The standard refining programs can work with twin 
models, but only if they consist of merohedral twins. In 



our case, the reflections do not correspond with an or- 
thorhombic unit cell, but with a, twice as large, pseudo 
cubic unit cell. Here we will explain the way we worked 
around the refining program. 

First, a hemisphere in hkl space was measured with 
-20 < h< 20, -20 < k < 20 and < I < 20. On this data 
set, we could refine our twin model with the six fractions, 
including the Ca concentration on the A site. Measuring 
such a large range in hkl space requires a large amount 
of time, roughly 10 days in the present set-up. Typically, 
for crystals with an orthorhombic unit cell measuring one 
octant of reciprocal space is sufficient. From refinements 
on selected parts of the data set, we concluded that we 
could investigate the structure, for instance at different 
temperatures, by measuring only the positive octant of 
the main fraction, thereby considerably limiting the mea- 
suring time. 

Here we describe the model and the application by 
SHELXL. The standard refining programs can work with 
merohedral twin models. However, in our case, the ob- 
served lattice does not coincide with one orthorhombic 
unit cell. Every observed reflection has contributions of 
up to six twin fractions. We measured both in the double 
cubic setting as in the orthorhombic setting of the main 
twin fraction. We transformed the h"k"l" indices to the 
six possible twin orientations. Three of the possibilities 
are given by 

h'k'l' = + 1"), k", ~(ft" - 1") (Bi) 
h'k'l' = 1(1" + k"), h", - k") (B2) 
h'k'l' = + h"), I", - h") (B3) 

The other three can be acquired by changing h' for V and 
V for -ft'. 

A new software program TWINSXL was developed to 
transform the standard data file, HKLS, by using the 
appropriate transformation matrices from a second in- 
put fild~3. The new data file is constituted of lines of 
ft, k, I, intensity plus standard deviation and twin fraction 
number. We used the "HKLF 5" option of SHELXL to 
refine data. The refinement uses a crystal model for the 
orthorhombic structure with the normal, adjustable vari- 
ables. Five variables for the twin fractions were added, 
the sixth fraction is calculated as the complement of the 
other five fractions. The sum of the appropriate calcu- 
lated intensities for all fractions was compared with the 
observed integrated intensities. 
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